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Rate-Based Model Predictive Control of Turbofan 

Engine Clearance 

Jonathan A. DeCastro 
QSS Group, Inc. 

Cleveland, Ohio 44135 

Abstract 

An innovative model predictive control strategy is developed for control of nonlinear 
aircraft propulsion systems and sub-systems. At the heart of the controller is a rate-based 
linear parameter-varying model that propagates the state derivatives across the prediction 
horizon, extending prediction fidelity to transient regimes where conventional models begin 
to lose validity. The new control law is applied to a demanding active clearance control 
application, where the objectives are to tightly regulate blade tip clearances and also 
anticipate and avoid detrimental blade-shroud rub occurrences by optimally maintaining a 
predefined minimum clearance. Simulation results verify that the rate-based controller is 
capable of satisfying the objectives during realistic flight scenarios where both a 
conventional Jacobian-based model predictive control law and an unconstrained linear- 
quadratic optimal controller are incapable of doing so. The controller is evaluated using a 
variety of different actuators, illustrating the efficacy and versatility of the control approach. 
It is concluded that the new strategy has promise for this and other nonlinear aerospace 
applications that place high importance on the attainment of control objectives during 
transient regimes. 


I. Introduction 

A LTHOUGH model predictive control (MPC) has been used for decades in some form or another in industrial 
processes, it is gaining ever more acceptance in aircraft propulsion applications due to advances in computing 
power of modern on-board control platforms. 1 2 Fast computation is necessary because the controller must generate 
current and future actuator commands based on open-loop receding horizon trajectories at each time instant. 3 The 
trajectory information allows the controller to optimize over the horizon while avoiding a pre-determined set of 
constraints, such as saturation limits of an actuator or operational limits that define the safety margins of a particular 
plant or process. In specific propulsion applications such as aircraft engine control or active turbine tip clearance 
control, this control strategy offers a wide range of utility. 

Figure 1 shows a general conceptual sketch of an active clearance control system, where the actuated shroud is 
depicted as responding to changes in clearance measured by a clearance probe. It is envisaged that replacement of 
open-loop clearance control systems aboard modern aircraft engines with faster actively-controlled systems will 
result in ultra tight clearances in the high-pressure turbine, permitting lower emissions and higher fuel savings by 
decreasing specific fuel consumption (SFC) as well as longer life through decreased exhaust gas temperatures 
(EGT). 4 Due to the wide variations in component deformations that occur during flight, if the clearance gap is not 
large enough, the turbine blades may mb against the shroud and permanently damage the blade tips. These 
variations arise due to the mismatch in thermal deformations and deformation rates of the rotor and stator assemblies 
along with centrifugal expansion and contraction of the rotor itself. However, on the basis that every 0.01 -inch 
reduction in clearance results in nearly 1% improvement in efficiency and 10°C reduction in peak EGT during 
transients, 5 the need to control clearance as tightly as possible becomes abundantly apparent. Certainly, advances in 
actuation and sensing hardware for active clearance control will allow for fast and accurate positioning of the shroud 
relative to the blades, 6 ' 7 but at these tight clearances, detrimental blade rubs may still arise due to the fact that there 
is no way of imposing mechanical limits on the clearance. The advantage of applying MPC to this application is 
that an inherent degree of robustness against blade mbs may enable operation with the tightest clearances that the 
hardware will allow, enabling the maximum operational benefits to be extracted. 
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Clearance reduction 
mechanisms: 

• rotor accel 

• rotor heating 

• case cooling 

• aero loads 


Clearance increase 
mechanisms: 

• rotor decel 

• rotor cooling 

• case heating 

• aero loads 

Figure 1. Active turbine clearance control concept (cutaway view of turbine blade tip and shroud). 

Because MPC is heavily reliant on an open-loop plant model, the model itself should be an accurate 
representation of the plant; in this case, a turbofan engine. In numerous cases, the pitfall in applying successful 
computationally-intense model-based control strategies is in the complexity of the model, as is the case with 
nonlinear models with a high degree of parameterization (e.g. neural nets). On the other hand, lack of fidelity can 
result in poor performance or robustness. In this paper, a gain-scheduled, linear parameter varying (LPV) MPC law 
is considered that retains both computational efficiency and prediction fidelity. In the past, gain-scheduling has 
often been applied under the assumption that analytical models of the plant exist and are at the disposal of the 
control engineer, as in Heise and Maeiejowski 8 and Shamma and Cloutier. 9 Unfortunately, because high-fidelity 
numerical engine models often used in gain scheduling applications are highly complex, it is nearly impossible to 
arrive at representative closed- form analytical models. A smaller group of papers, such as Mehra et a/, 10 treat the 
use of Jacobian linearized LPV models for MPC, but the limitation of traditional LPV models is in the loss in 
prediction quality when transients cause the models to move away from quasi-equilibrium. Over the horizon, these 
models may substantially depart from the near-equilibrium regions of the parameter space where model validity 
holds and substantially compromise optimality. Leveraged by a rate-based linearization method introduced by Leith 
and Leithead," an alternate framework conducive to rapidly-varying plants is explored in this work, which operates 
upon the state derivatives instead of the states themselves. A virtue of this transformation is that fidelity during 
transients is preserved, enhancing the prediction and global performance of the MPC. To the author’s knowledge, 
this is the first attempt at developing a rate -based LPV modeling technique (hereafter abbreviated as R-LPV) for 
MPC. The results of the assessment in this paper are useful as further guidelines from which to build MPC 
controllers for other aerospace applications (e.g. engine control), in addition to the application described here. 

This paper is organized as follows. To lay the groundwork for the new MPC, a description of both conventional 
and rate-based linearization methods for LPV representations as applied to a numerical turbofan engine model are 
presented. Next, a description of an LPV model predictive controller is given, along with the proposed disturbance 
model estimation structures. Aspects of applying the MPC to a fast active clearance control device are then 
described, and details of the actuators used and the high fidelity engine model are given. Finally, results of active 
clearance control simulations are presented, demonstrating the comparative advantages of the proposed controller 
over the conventional approach. In that section, both servohydraulic and thermal deformation-based actuators are 
investigated to illustrate the versatility and robustness of this control approach. 

II. LPV Models of Nonlinear Processes 

The purpose of this section is to introduce concepts relating to linearization of numerical models in the LPV 
framework, as motivated by the highly nonlinear nature of the turbofan engine. Once these methods are established, 
germane properties of the two methods are described, as they apply to model predictive control. 

General nonlinear systems of the form x = f(x,u,v ) , y = g(x,u, v) can be placed in the following standard LPV 
form: 

x = A(p)x + B t {p V + B 2 {p)v 1 

y = c(p)x + o, (p)i + d 2 (p)v 


clearance gap 


Active clearance 
control 
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The state vector, control input, disturbance input, and output vector are xeR" , u e R "“ , v e R" , and yeR" y , 
respectively. The scheduling parameter vector is denoted by p and can be either exogenous, an input external to 
the plant, or endogenous, a state of the system. In the latter case, the system is said to be quasi-LPV. 

Two Perspectives on Linearization 

Jacobian linearization 12 is often considered the default method for attaining representations of nonlinear 
processes, especially those for which analytical expressions do not exist. Although it is possible in theory to 
linearize away from equilibrium, practical considerations often limit perturbation of the model or actual system to 
steady-state operating points. This procedure, when conducted at several distinct operating points, produces a 
family of linear models 


& = A(p)Sx + B t (p\ + B 2 {p\ 5v 
Sy = C(p)Sx + Z), ( p)Sa +D 2 {p)Sv 

where the perturbation values 5x = x - x 0 , Sii = u -u 0 , Sv = v - v 0 , and Sy = y - y 0 are functions of the steady-state 
trim values, x 0 , u 0 , v 0 , and y 0 . Although this family can be quite large when adequate coverage of the entire 
operating envelope is needed, the result is by no means global as the state derivative trim value x 0 must still be 
assumed to remain zero everywhere. This implies that LPV trajectories are limited to a subset of actual trajectories 
of the nonlinear system exactly at equilibrium, but this restriction is often relaxed by permitting transitions between 
points so long as variations in p are sufficiently slow. Within the neighborhood of a single linearization point, the 
LPV input/output relationship is accurate to within first-order of the actual system. 

Practically speaking, the family of state-space models is scheduled by either a nearest neighbor or interpolation 
approach. Since the LPV system in Eq. (2) is actually a nonlinear affine system, trim value updates must be 
included along with the system matrix updates, as depicted in Fig. 2. If the nonlinearity is particularly dominant, i.e. 
highly sensitive to minor variations in p , then quasi-LPV representations may be rendered invalid even with slow 
variations due to the fact that the trim values are only approximations as the trajectory progresses between operating 
points. Moreover, the recursive nature of the quasi-LPV parameter update causes any error to accumulate without 
bound, possibly leading to significant model/plant disparities even within small time frames. 

In an alternate linearization perspective known as rate-based linearization, 11 the nonlinear system is 
approximated to first order by an affine Taylor series expansion as before, but the states are transformed by simply 
taking their derivatives. In so doing, the assumed static first-order terms disappear, leaving the actual system state 
derivatives as the LPV states. Perturbation values now become unnecessary, and the system of equations takes on 
the following form: 


x = A(p)x + B t (p)u + B 2 (p)v 
y = C{p)x + £>, {p)u + D 2 (p)v 


In this representation, the state vector is now represented by £, = [i 7 y 7 ] and the state matrices are augmented 
appropriately. To complete the R-LPV model, initial states must be specified if the model is to start at points other 
than the origin. Figure 3 shows an input/output signal diagram for this type of model. If the trim values are, in fact, 
obtained at steady-state, it may be observed that the perturbation states Sx in Eq. (2) are first-order approximations 


«o ip) 
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Figure 2. Jacobian LPV models. 
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Figure 3. R-LPV models. 

to the derivatives of the actual states x in the limit sense. The family of Jacobians obtained in Eq. (2) may thus be 
substituted into Eq. (3) and scheduled identically. In contrast to Eq. (2), however, all trim values including the state 
derivative trim x 0 become implicit to the model. To elaborate, if p is sufficiently slow during a particular 
transient, then the nonhomogeneous terms in Eq. (2), a linear combination of trim values obtained at equilibrium 
that are dependent upon p , will be close enough to reality that the resulting trajectory will closely correspond to the 
actual. If the system were excited by another transient where p is now fast, the nonhomogeneous terms will be 
incorrectly computed from the subset of trim values at equilibrium, causing inconsistencies between the model and 
actual trajectories. In the R-LPV setting, these nonhomogeneous terms are assumed constant, resulting in Eq. (3) 
being linear. No restrictions are therefore imposed on these terms being constructed from equilibrium trim values. 

R-LPV models have improved fidelity at all neighboring points in the parameter space, including non- 
equilibrium conditions, implying that the slow variation condition may be relaxed to some extent. Furthermore, by 
removing the state derivative trim and all remaining trim values from the formulation, these quantities no longer 
need to be approximated. This does not solve the problem of obtaining new Jacobians away from equilibrium, but 
the reduced parameterization and interpolation thereof lessens propagation of error across the trajectory existent in 
quasi-LPV formulations. In the context of receding horizon control, improved fidelity means that global optimality 
is more closely preserved. Because the open-loop model is effectively a quasi-LPV model that recursively executes 
several samples into the future, reduced error propagation becomes equally important in realizing global optimality 
as maintaining fidelity in non-equilibrium regimes. 

The penalty for the improved model is that the state vector dimension is now expanded to n + n y and that the 
computed control signal is the integral of the actual controller output, u = \u ' , due to the rate-dependency of the 
controller structure. Conversely, the measured disturbance inputs to the LPV models are brought in as a derivative, 
i.e. v' = dv/dt . It should be cautioned that introduction of uncertainties may certainly have a deleterious effect due 
to this derivative expression if these inputs are susceptible to noise or biasing, especially in plants with high input- 
output sensitivities. In practice, this disturbance is often manifested as a demand-related event that is subject to little 
noise or other types of uncertainty. Moreover, estimation structures that operate upon separate disturbance states 
decouple the estimation loop from these effects. 

III. Model Predictive Controller in Rate-Based LPV Form 

In this section, the practical problem setup for MPC is given. First, the disturbance model structure for output 
estimation in rate -based form is presented with treatment given to consolidating the estimator states with the 
expanded R-LPV states. Then, with this complete plant/estimator model, the R-LPV formulation of MPC is 
provided. 

Disturbance Models in the Rate-Based Paradigm 

In receding horizon control, input/output disturbance models are necessary to reduce the effects of any plant- 
model mismatch and to absorb exogenous disturbances affecting the plant but unknown to the model. Due to lack of 
information regarding most disturbances, such models are often chosen to be pure integrators so that any constant 
disturbance may be asymptotically rejected and output tracking may theoretically occur without error. It is well 
known that, if the augmented system is detectable and the number of disturbance states equals the number of 
outputs, then all output channels will have zero steady-state error. 1314 Normally, the state vector is expanded to 
incorporate these new disturbance states, d e R" . However, this approach can be improved to avoid undue 
computational complexity arising from incorporating an additional n y integrator states (a total of 2 n v more than the 
original model). Recognizing that the output integration is decoupled from the plant dynamics, the disturbance 
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model inputs are simply added to y in Fig. 3 and the state sub-vector y in Eq. (3) now becomes y + d . The state 
space equations are then 


A = A ^ + B { u' + B 2 v' + Bpv p 
y = C^ + w m 


( 4 ) 


with the augmented matrices 


A = 


4p) 

c(p) 


B t = 


bAp) 

dAp) 


b 2 = 


B i{p) 

d 2 (p) 


B,= 


C = 0 I , 


As in Eq. (3), the augmented state vector E, is composed of the state derivative vector and the output vector. 
Uncorrelated additive process and measurement noise intensities, xv p e R" and w m g R" , enter the disturbance 
model input channels and measurement channels, respectively. The diagonal fdter matrix A n is included to avoid 
pure integration at the outputs. Note that the R-LPV realization augmented with disturbance models will have the 
same number of augmented states as conventional LPV realizations with disturbance models. 

Model Predictive Control Law 

Given the R-LPV plant/estimator model, optimal control actions are found by minimization of a finite-horizon 
linear quadratic objective function with respect to the controls and the outputs. 3 Converting to discrete time, the 
linear-quadratic cost function for the MPC law may be written as 


J k = X IK-+.- -yrMifoifk+i -y r Mi) +u Li Ru k +i . 


( 5 ) 


subject to the discretized version of the plant model of Eq. (4) and constraints on the manipulated controls and 
measured outputs at time k. The subscript r denotes the reference setting and p represents the prediction horizon. 
The symmetric, positive-definite weight matrices Q and R are to-be-determined design parameters, and the usual 
requirement holds that the pairs ( A,B t ) and [a,Q 1 2 ) are stabilizable. A Kalman filter state estimator is designed 
on the condition that (c,a) is detectable, assuming the process and measurement noise covariances W and V, 
defined by the intensities w p and vv„„ respectively. It is important to note that, in order to impose any weighting or 
constraints on the actual controls in the rate -based MPC framework, the output vector must be augmented with those 
control inputs. 

The LPV MPC formulation is achieved by computing state prediction equations recursively at each time step. 
Using the same state space notation from continuous time for the discrete time system equations, the state 
predictions are computed by 


4k + 1 


1 

<1 

i 


V* 


= A£+BX_ 1+ Bf 


+ B 2 


1 

+ 

"a 

1 


> 

+ 

1 


y'k +P - 1_ 


(6) 


where E, k are the estimated R-LPV states. The prediction matrices are 


Aik ) 


BAk) 


, B 1 = 


IF(*+0 

_ i=0 


YjFAj,k + p)+BAk + p-l) 

. J =1 
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Bf = 


B\k) 


Y j F' l {j,k + p)+B { {k + p-\) Y j F x (j,k + p)+B l {k + p-\) ■■■ B t {k + p- 1) 

M 


P-2 

I 

7=1 


B, 


S 2 (/t) 0 0 

77, (jo -!,£ + /?) F 2 (p-2,k + p) ••• F 2 (l,£ + p) 


where 








In this representation, the control inputs are split into the present value of the control derivative u’ k and a vector of 
tut li re control derivative moves. 

The trajectory of predicted outputs is 


y k+ i =c(k + i% +i +D 2 (k + i)v' k+i (7) 

for each i = l,...,p . The dependence of the Jacobians on p is implied by their discrete-time dependence on k. 
From the above equations, it is readily apparent that the internal models assume quasi-LPV formulations (an 
internally-computed scheduler) in forming the predictions, but may use an external scheduler for state estimation as 
long as state continuity between each realization is guaranteed. Because the internal model must be strictly proper; 
i.e. zero gain at infinite frequency, it is necessary to formulate Jacobian LPV systems without the D x term. In 
contrast, the fundamental model in R-LPV systems may be of any type because the R-LPV transformation is, by 
definition, one that maps any arbitrary system to a strictly proper system because the original D matrices become 
part of the B matrices. 

Standard quadratic programming (QP) approaches are available to efficiently solve this convex constrained 
optimization problem; the approach employed here uses an “active set” method that optimizes on a subset of 
inequality constraints that are active at the current step. Common ways of making the control law computation more 
efficient are to limit the control horizon m to some m < p or to update only intermittently (known as blocking). To 
expand the feasible set of the QP problem, the constraints are not posed as hard limits but instead are softened by 
incorporating a slack variable on the cost function that introduces an expensive penalty should any constraint 
violations occur. Computation is further streamlined by generating the Jacobians, estimator gains, and prediction 
matrices in an offline batch mode as a set of parsed arrays so that the on-line algorithm only needs to perform linear 
interpolation on these arrays at each time step. 

IV. Active Tip Clearance Control in a Commercial Turbofan Engine 

As depicted in Fig. 1, the clearance control actuator is designed to move a segmented shroud structure toward 
and away from the turbine blades. The concept, as referred to in Lattime, et a/ 15 consists of several identical 
actuators that are placed at locations around the shroud circumference with each effecting movement in one shroud 
segment. It is assumed that each actuation point has a collocated clearance measurement at the shroud’s frame of 
reference. In this configuration, the detected clearance is the difference between the combined deformation of 
turbine components (referred to as the un-actuated clearance) and the actuator displacement, where a positive 
displacement results in smaller clearance. The active clearance control system can operate either upon a minimum 
detected clearance and single command signal, or upon several local clearance signals and multiple actuator 
commands. The MPC is adaptable to both configurations, but only the single-output case is considered here. 
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Two actuator platforms were used to demonstrate predictive control behavior across different realistic hardware 
design scenarios: one representing a servohydraulic actuator and the other representing a thermal deformation-based 
actuator. A detailed nonlinear model of the servohydraulic actuator described in DeCastro and Melcher 16 is 
implemented as the truth model and a parsimonious reduced order model is used in the MPC. The dynamics of the 
actuator are primarily driven by the servo control loop and are approximated here by the following third-order 
system with corner frequency of approximately 0.2 Hz: 


1283 

(s + 222>2\s + 4.035X-S- + 1.424) 


( 8 ) 


The actuator stroke is measured by a transducer for servo-control and is therefore accessible to the engine control 
unit as a measured state of the system. The second actuator, a thermal-deformation based device, is similar to 
current thermal actuation systems that utilize secondary engine air to regulate the flange attachment to the shroud 
structure by way of a proportional air valve. 4 This device is represented by the first order system 



zs + 1 


(9) 


Because the intent of this case is purely illustrative, the thermal actuator internal model and truth model are taken to 
be identical. For both actuator types, constraints were placed on the minimum and maximum allowable control 
input settings, corresponding to zero and 0.1 -inch displacements, respectively. 

For this tip clearance application, a commercial engine simulation (CES) representative of a large, high-bypass 
turbofan engine coupled with the actuator model was used as the baseline truth model for closed-loop evaluation. 
The simulation features a detailed model of the turbine section component deflection to accurately represent 
demand-related clearance effects. On the basis that the MPC is to be implemented as an isolated control loop from 
the engine controller, a simplified closed-loop engine fan speed compensator was incorporated. The compensator is 
a linear, second order system with a free integrator at the output, designed to obtain a slightly overdamped nominal 
closed-loop bandwidth of approximately 0.6 Hz. The requested fan speed is actually a mapping function consisting 
of power lever angle ( PLA ) and three environmental parameters (altitude, Mach number, ambient temperature). 
Multi-mode controllers that typically deal with limits normally encountered with large transients have not been 
incorporated, but the MPC structure can be viably altered in practice to accommodate this additional level of 
sophistication. Alternately, if the tip clearance MPC is implemented in a unified engine MPC, separate models of 
the controller are unnecessary because the MPC generates these engine commands. 

Nine linear models were generated of the closed-loop engine at equispaced PLA s between ground idle and 
maximum power at sea-level static conditions, constructed with several sensed output channels and two input 
channels: delta clearance (provided by the actuator) and fuel flow. Jacobians at these points were found by 
perturbing these inputs, nine engine states (fan and core shaft velocities and seven internal metal temperatures), and 
the three fan speed controller states. The turbine component states are numerous and were omitted from the 
linearization. Because the operating conditions considered here are limited to constant altitude and airspeed, the 
scheduler is taken simply as the measured fan speed, XN2. 

With the hydraulic actuator model, the LPV model has a total of two inputs, fifteen states, and six outputs; these 
I/O’s are listed in Table 1. The total number of R-LPV internal model states is 22 (n + n v +«„), only one state 
greater than the Jacobian LPV counterpart. Figure 4 shows the interconnection of the rate-based active clearance 
control MPC and the turbofan engine plant. The exogenous disturbance ( 8 k ) is brought in only to the engine, while 
the demand-induced measured disturbance ( v k ) input, PLA, is brought in to the engine and MPC. This may be 
extended to capture other sensed inputs (ambient temperature, Mach number, altitude) if needed. As mentioned 
earlier, differentiation of input variables is necessary in the rate-based internal model framework. Such a 
transformation can amplify any noise or biasing on the input channels, but since this application calls for a 
disturbance that is a function of throttle position, it may be argued that the signal is mostly free of noise and other 
sources of uncertainty. Nonetheless, examining the effect that these uncertainties have on transient performance is 
warranted to ensure that optimality is not compromised. 
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Table 1. Inputs and outputs for R-LPV tip clearance controller. 


Inputs 

Outputs 

Actuator command, UACT 
Power lever angle, PLA 

Fan speed, XN2 
Fuel flow, WF 
Exhaust gas temperature, T5 
Stage 1 clearance, CLR 
Actuator position, XACT 
Actuator command, UACT 


Design Considerations 

In this tip clearance implementation, minimization of SFC and EGT is indirectly achieved by penalizing 
clearance excursions away from a small-valued set point. Hence, the objective function and loop transfer function 
reduce to single-input, single-output (SISO), with the actuator command as the control input and clearance as the 
regulated output. The unconstrained SISO control input is then used exclusively for design: 

u[=-Kl (10) 


where K is the full-state feedback gain matrix. 

The tuning procedure selected in this paper follows Doyle and Stein, 17 where the objective is to first specify a 
desired loop transfer function (closed or open), then apply a quasi-deadbeat estimator to the loop to asymptotically 
achieve the specification. With the ratio R/Q set to 0.001, the controller attains fast response but avoids high- 
frequency gain/phase distortions in the control loop return ratio K (zl - A ) 5, due to discretization effects. In so 
doing, the return ratio more closely resembles its continuous-time counterpart, which has the property of providing 
infinite gain margin against gain increases and 60° of phase margin. 

V. Demonstrative Simulations 

In this section, simulated results of the active clearance control MPC application are given and comparisons 
made between controllers derived using conventional internal modeling methods and the rate-based method 
explained in this paper. Additionally, ties are made to the practicality of this MPC implementation in active 
clearance control using the two actuation methods described earlier. For all reported results, a 50-Hz sample rate 
was used. The minimum clearance constraint and clearance set point were set to identical values, which presents the 


MPC 



Figure 4. Block diagram of MPC with R-LPV models. 
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most demanding scenario for the controller since maintenance of the clearance set point and avoiding blade rubs are 
now invoked as competing objectives. 

Distinctions between Rate-Based and Jacobian MPC 

Clearance responses with the Jacobian-MPC law for the servohydraulic actuator given a large-magnitude PLA 
transient are shown in Fig. 5. For this assessment, the prediction horizon p was set to 20 samples and the control 
horizon m to 2 samples, with both these settings chosen to obtain an unconstrained response that resembles the 
infinite-horizon response. In the figure, fan speed and adjusted clearance (top and bottom traces, respectively) 
represent measured outputs. While not a physical output, the un-actuated clearance shown in the middle trace 
represents the combined deformation of the turbine section components, a demand-related effect from the engine 
essential for predicting the evolution of clearance through time. This quantity is computed as the sum of the 
controlled clearance and the actuator displacement (not shown), which is a valid estimate only since the feedback 
effects of clearance on the engine’s gas path dynamics are considered to be minor. 

In the bottom trace of Fig. 5, it is clear that the transient causes the adjusted clearance to fall below the constraint 
(dashed line) at about 1.8 seconds, which is in violation of the objective that the adjusted clearance must be kept 
positive to avoid the possibility of blade rubs. The open-loop prediction trajectory overlays reveal that the model 
incorrectly predicts the evolution of fan speed and un-actuated clearance outputs, especially during the time interval 
between 1.8 and 2.2 seconds, which correlates to the time when the constraint violation begins. The output 
trajectory plot in Fig. 6 is helpful for visualizing gradient information in the output space, containing the same 
prediction trajectories as Fig. 5 plus one additional trajectory at the 1.8-second mark. As expected, the 1.8-second 
prediction has the largest initial gradient error of the four trajectories, as shown by the large migration away from 
the actual trajectory early on in the horizon. 

The improvements that the R-LPV models afford the optimization are shown in Fig. 7, given the same PLA 
transient and an identically-tuned controller (i.e., one with the same design parameters). In comparison with the 
previous case, the adjusted clearance (bottom trace) stays at or above the constraint, satisfying the blade rub 
objective. The prediction traces in Fig. 7 and associated output trajectory plot in Fig. 8 show that the sharp, erratic 
responses disappear, with initial gradients following much closer to the actual. Note that the 1.5-second prediction 
traces settle out at constant rates-of-change. The reason for this is because the measured disturbance derivative, not 
the measured disturbance itself, is held constant throughout the horizon in the R-LPV setting. For all simulations 






Figure 5. Jacobian LPV responses, showing actual 
(red) and predicted outputs (black dots). 


Figure 6. Jacobian LPV output trajectories, 
showing actual (red) and predicted outputs (black 
dots). 
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Figure 7. R-LPV responses, showing actual (red) 
and predicted outputs (black dots). 


Figure 8. R-LPV output trajectories, showing 
actual (red) and predicted outputs (black dots). 


executed, this does not negatively impact the convergence of the controller. At 3.0 seconds, all internal model 
transients die out, with any subsequent error being a result of the slow neglected clearance dynamics. The closer 
agreement between the R-LPV models and the nonlinear plant throughout the transient is owed to the fact that the 
non-homogeneous terms in the R-LPV models are no longer explicitly parameterized and therefore not subject to the 
associated limitations. 


Practical Assessment of the Rate-Based MPC 

Where the previous evaluation was only concerned with utilizing the feedforward attributes of the MPC in 
optimizing control actions based on demand-induced clearance events, actual in-flight situations require the need to 
reject exogenous disturbances and maintain performance with changing parameters, an exercise that is largely a 
question of feedback compensation. In the active clearance control application, clearance changes induced by flight 
loads acting on the engine are considered, whose effects become more problematic as clearance set points decrease 
because of the decreased margins against blade mbs. As alluded to in Olsson and Martin, 18 there are multitudes of 
engine transients that occur during flight that are of interest in active clearance control; a sampling of which includes 
takeoff, re-burst, thrust reversal, hard turns, and aircraft stall. Each event contributes a demand-induced component 
on clearance due to PL A, etc., and a presumed unknown flight load-induced component as a result of turbulence, 
maneuvering of the aircraft, and asymmetric thrust loads that cause the engine to bend. Because of their large, 
rapidly-acting contributions on clearance, the takeoff, thrust reversal, and aircraft stall events are employed to 
exercise the MPC laws given realistic clearance events that occur during normal operation; the magnitudes and 
rates-of-change of each event were reproduced approximately from actual flight data reported in Ref 18. Since 
takeoff and re-burst are essentially identical in nature aside from the thermal state of the turbine at the start of the 
event, examining the transient MPC response with a takeoff event effectively serves the dual purpose. 

To compare the MPC strategy with a controller very similar to a linear quadratic Gaussian (LQG) control law, 
the clearance constraint was removed from the MPC computations, rendering the new controller to be similar by 
nature to a gain-scheduled LQ controller with anti-windup protection. The purpose of this comparison is to justify 
the benefit of MPC in the face of the additional computational burden. Note that the unconstrained, LQG-like 
controller is still globally valid in the sense of Eq. (3). If a truly gain-scheduled LQG controller were used, this 
would be optimal only in the vicinity of the linearization points. Figure 9 shows the MPC response given the three 
event scenarios at nominal conditions for a thermal actuator with 1.6-second time constant (r = 1.6). This case 
represents an actuator whose response rate roughly matches the maximum deformation rate of the rotor during a 
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Figure 9. R-LPV responses for fast thermal actuator. 

power-increase event. To extend the anticipation capability, the prediction horizon p was increased to 50 samples 
and control horizon m to 5 samples. From the figure, the MPC response to a thrust reversal transient produces the 
most benefit over an unconstrained controller, eliminating the 2-mil excursion below the constraint. It is 
encouraging that enforcement of the clearance constraint is a result of very smooth actuator command signals rather 
than abrupt or highly oscillatory control actions. Although unconstrained responses in the remaining two events 
have milder constraint violations, it is important to note that the constrained controller satisfies the blade rub 
objective in both cases. One may interpret these results from a systems point-of-view. If the clearance setpoint 
using this LQG-like controller must be 2 mils greater than the MPC setpoint to accommodate the possibility of the 
worst-case event, then this would detract from exploiting the benefits of a tip clearance control system. Based upon 
past studies, 5 the penalty on efficiency would be 0.2% and the penalty on EGT overshoot would be 2°C. 

Figure 10 shows the response with a thermal actuator with r = 5.3 , a rate 60% slower than the previous case. 
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Figure 10. R-LPV responses for slow thermal actuator. 
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Table 2. Minimum clearance for constrained and unconstrained R-LPV MPC. Values are in mils 



(1x10- 

3 in). 



Unconstrained Constrained 

Delta 

Thermal Actuator, r = 

1.6, nominal engine 


Takeoff 

-0.28 

-0.01 

0.27 

Thrust reverse 

-1.93 

-0.01 

1.92 

Aircraft stall 

-1.00 

-0.04 

0.96 

Thermal Actuator, r = 

5.3, nominal engine 


Takeoff 

-5.01 

-4.69 

0.32 

Thrust reverse 

-10.5 

-9.80 

0.70 

Aircraft stall 

-8.71 

-8.02 

0.69 


Hydraulic Actuator, nominal engine 


Takeoff 

-0.80 

-0.18 

0.62 

Thrust reverse 

-2.29 

-0.50 

1.79 

Aircraft stall 

-3.26 

-0.71 

2.55 


Hydraulic Actuator, off-nominal engine 


Takeoff 

-0.60 

-0.13 

0.47 

Thrust reverse 

-2.15 

-0.50 

1.65 

Aircraft stall 

-3.26 

-0.70 

2.56 


This case is included to illustrate the consequences of applying MPC to an actuator that is slower than the rate of 
mechanical rotor growth, similar to state-of-the-art actuators. In the figure, the range of the control signal lies 
between one and zero and saturates beyond those values. For both controllers, the control signal saturates at zero in 
attempt to move the shroud away from the blades. With the constrained controller, anticipation of the event in the 
takeoff and thrust reverse events results in a slight increase in adjusted clearance before becoming altogether 
negative upon loss of authority. Although no anticipation is possible in the purely exogenous stall event, the 
adjusted clearance stays at zero longer than for the unconstrained case, a result of the more rapid compensation of 
the constrained MPC when riding the clearance constraint. It may be suggested that MPC performance can be 
further improved with anticipation of the exogenous events, perhaps either by direct estimation of the flight load 
event or by correlating to other measurements (e.g. structural sensors). These effects may not be assumed uniform 
around the turbine circumference. 

Table 2 summarizes the MPC benefits for the three actuators and also presents a case of tip clearance control 
amidst off-nominal engine conditions. Note that any violations of the constraint are a result of the slack tuning 
parameter described earlier. Comparing the fast thermal actuator with the nominal hydraulic device, both the 
unconstrained and constrained clearance values are quite different, which is directly related to sub-optimality 
introduced with plant/model mismatch. Recall that the nonlinear hydraulic actuator is adopted as the truth model, 
while the thermal actuator only uses the linear model. With the exception of the slow thermal device, enforcement 
of the constraint results in commensurate amounts of clearance improvements for each. As echoed in Fig. 10, the 
MPC applied to the slow actuator offers some improvement, albeit less than a full-authority actuator does. The 
degraded engine case is implemented as a shift in engine parameters to roughly represent 6000 cycles of engine 
service life degradation, as described in Sallee, 19 coupled with a shift in fuel flow sensitivity of -5%. In a gross 
dynamical sense, any engine degradation is not expected to significantly impact the behavior of the engine, and 
much less the behavior of clearance dynamics. This important attribute is exemplified by the results. 

To address the aforementioned issue concerning the integrity of the internal model amidst uncertainties in the 
measured disturbance signal, simulation cases were generated where the PLA measured disturbance input was fouled 
by uncertainties. The results are shown in Fig. 11, where two cases show the effects of the MPC receiving PLA 
commands corrupted by a ±5% shift in sensitivity and a third where the PLA has additive white noise with variance 
of approximately 6% of full scale. The traces show that the introduction of noise causes excitation of higher- 
frequency controller modes, but no adverse effects in control performance. Although the controller sees a PLA 
derivative that differs greatly from the actual value entering the engine, the deadbeat estimator is quick to reject 
these disparities and the prediction is largely unaffected during the comparatively short span of the horizon. 
Naturally, longer horizons will exacerbate the effects of noise on the input channel, possibly producing an erratic or 
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Figure 11. Constrained R-LPV responses for a thrust reverse transient with corrupted measured 

disturbance. 

even divergent control signal. Nevertheless, these results are important in confirming that R-LPV models are indeed 
viable for use in MPC even in the presence of noise. 

VI. Conclusions 

In this paper, an innovative model predictive control (MPC) strategy using rate-based linear parameter varying 
(R-LPV) models is introduced. The virtues of the rate-based framework permit state derivatives to be directly 
specified, resulting in the implicit inclusion of state derivative trim values in the family of trajectories within the 
neighborhood of linearization equilibrium points. This enhances transient performance by widening the model’s 
valid applicable range without introducing additional linearization points, especially at non-equilibrium. Removal 
of the explicit trim value information furthermore results in models that exhibit less approximation error, an 
important attribute of any recursive modeling technique, such as quasi-LPV. 

The new MPC strategy is applied to an active clearance control system, where it is highly important to minimize 
clearance but avoid violating the minimum clearance constraint that protects against detrimental turbine blade mbs. 
In this work, the active clearance control optimization problem consists only of maintaining a minimal clearance set 
point, however directly optimizing against performance metrics such as exhaust gas temperature ( T5 ) and fuel flow 
(WF) is a straightforward extension using this framework. The R-LPV strategy was demonstrated to maintain 
prediction fidelity and control optimality in transient regimes where the conventional Jacobian LPV method does 
not. Additionally, the controller was proven to be advantageous for tip clearance control when comparing against an 
analogous unconstrained, gain-scheduled controller similar to LQG. 

With regard to computational complexity, although the R-LPV transformation produces models with state 
dimensions greater than conventional LPV by the number of controlled inputs and outputs, it is shown that the 
integrated output states may be merged with the disturbance model states. As such, the internal model states are 
greater by only the number of inputs, which is usually small in propulsion system applications. Although the new 
MPC strategy was applied to the rather demanding active clearance control problem with competing objectives, the 
theoretical groundwork presented here provides sufficient generality to carry over to MPC applications in other 
nonlinear control problems; among the possibilities being turbofan engine control. Because of the demonstrated 
advantages this method has over conventional methods, it has great potential for control of aerospace applications 
where maintaining control optimality during transients is especially critical. 
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